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The masses of heavy tetraquarks with hidden charm and bottom are calculated in 
the framework of the relativistic quark model. The tetraquark is considered as the 
bound state of a heavy-light diquark and antidiquark. The light quark in a heavy- 
light diquark is treated completely relativistically. The internal structure of the 
diquark is taken into account by calculating the diquark-gluon form factor in terms 
of the diquark wave functions. New experimental data on charmonium-like states 
above open charm threshold are discussed. The obtained results indicate that the 
X(3872) can be the tetraquark state with hidden charm. The masses of ground state 
tetraquarks with hidden bottom are found to be below the open bottom threshold. 

PACS numbers: 12.40.Yx, 14.40.Gx, 12.39.Ki 


Recently, significant progress in experimental investigations of charmonium spectroscopy 
has been achieved. Several new states A"(3872), Z(3931), E(3943), X(3943), E(4260) were 
observed jlj which provide challenge to the theory, since not all of them can be easily 
accommodated as the cc-charmonium. The most natural conjecture is that these states are 
the multiquark composite systems considered long ago e.g. in fl]. The proposal to revisit 
the multiquark picture using diquarks has been raised by Jaffe and Wilczek |3]. Currently 
the best established state is the narrow A"(3872) which was originally discovered in B decays 
and later confirmed in pp collisions j|6, jj. Its mass and observed decays, which favour 


J PC = 1 ++ assignment, make a cc interpretation problematic |8|. Different theoretical 
interpretations of the A"(3872) state were put forward which use the near proximity of its 
mass to the D°D*° threshold. The most popular ones are: the D° — D*° molecular state 
bound by pion and quark exchanges an s-wave cusp at D°D*° threshold Q and the 
diquark-antidiquark [ cq ] [cq] tetraquark state Q (q = u, cl). 


Maiani et al. 0 hi the framework of the phenomenological constituent quark model 


considered the masses of hidden charm diquark-antidiquark states in terms of the constituent 
diquark mass and spin-spin interactions. They identified the A(3872) with the S'-wave bound 
state of a spin one and spin zero diquark and antidiquark with the symmetric cliquark- 
spin distribution ([cg]s = i[cg]s = o + [cg]g =0 [cg]s = i) and used its mass to fix the constituent 
diquark mass. Spin-spin couplings were fixed from the analysis of the observed meson and 
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baryon masses. On this basis they predicted the existence of a 2 ++ state [cq] s i =i[cg],s=i that 
can be associated to the F(3943). They also argued [12] that F(4260) could be the first 
orbital excitation of the charm-strange diquark-antidiquark state ([cs]s =0 [cs]s =0 )p-wave- In 
Ref. 0 it is pointed out that non-leptonic B decays provide a favourable environment for 
the production of hidden charm diquark-antidiquark bound states. In contrast it is argued 
0 that the observed X(3872) production in B decays and in high-energy pp collisions is 
too large for a loosely bound molecule (with binding energy of 1 MeV or less). 

In this paper we use the relativistic quark model jEj |l(J based on the quasipotential 
approach to calculate the mass spectra of tetraquarks with hidden charm and bottom as the 
heavy-light diquark-antidiquark bound states ([Qq][Qq], Q = c,b). Recently we considered 
the mass spectra of doubly heavy ( QQq ) 0 and heavy ( qqQ ) [18] baryons in the heavy- 
diquark-light-quark and light-diquark-heavy-quark approximations, respectively. The light 
quarks and light diquarks were treated completely relativistically. The internal structure 
of the light and heavy diquarks was taken into account by calculating diquark-gluon form 
factors on the basis of the determined diquark wave functions. The found good agreement 
[18] with available experimental data gives additional motivation for considering diquarks 
as reasonable building blocks of hadrons. It is important to note that all parameters of our 
model were determined from the previous considerations of meson mass spectra and decays, 
and we will keep them fixed in the following analysis of heavy tetraquarks. 

In the quasipotential approach and diquark-antidiquark picture of heavy tetraquarks the 
interaction of two quarks in a diquark and the heavy diquark-antidiquark interaction in 
a tetraquark are described by the diquark wave function (T rf ) of the bound quark-quark 
state and by the tetraquark wave function (Tp) of the bound diquark-antidiquark state 
respectively, which satisfy the quasipotential equation of the Schrodinger type 


/ b 2 (M) p 2 ' 

\ 2 p R 2 p R 


d'd.'Hp) — 


d 3 q 

w? 


I/(p,q; q), 


( 1 ) 


where the relativistic reduced mass is 
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and Ei, E 2 are given by 


M 2 — ml + m\ M 2 — m\ + ml 

E\ — - 7T77 -> E 2 = — 


2 M 


2 M 


( 3 ) 


here M — E x + E 2 is the bound state mass (diquark or tetraquark), m^ 2 are the masses of 
quarks (q x and q 2 ) which form the diquark or of the diquark ( d ) and antiquark (d') which 
form the heavy tetraquark (T), and p is their relative momentum. I 11 the center of mass 
system the relative momentum squared on mass shell reads 

h 2 (M) _ W 2 ~ (mi + m 2 ) 2 ][M 2 - (mi - m 2 ) 2 } 
y J 4 M 2 ' 1 ; 

The kernel V(p, q;M) in Eq. (JTJ is the quasipotential operator of the quark-quark or 
diquark-antidiquark interaction. It is constructed with the help of the off-mass-shcll scatter¬ 
ing amplitude, projected onto the positive energy states. In the following analysis we closely 
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follow the similar construction of the quark-antiquark interaction in mesons which were ex¬ 
tensively studied in our relativistic quark model [jjJ[l!§]. For the quark-quark interaction 
in a diquark we use the relation V qq = V qq /2 arising under the assumption about the octet 
structure of the interaction from the difference in the qq and qq colour states . 1 An important 
role in this construction is played by the Lorentz-structure of the confining interaction. In 
our analysis of mesons while constructing the quasipotential of the quark-antiquark interac¬ 
tion, we adopted that the effective interaction is the sum of the usual one-gluon exchange 
term with the mixture of long-range vector and scalar linear confining potentials, where 
the vector confining potential contains the Pauli terms. We use the same conventions for 
the construction of the quark-quark and diquark-antidiquark interactions in the tetraquark. 
The quasipotential is then defined as follows [16, [l7j 
(a) for the quark-quark (Qq) interaction 

V (p, q; M) = Ui(p)u 2 (-p)V(p, q; M)m(q)u 2 (-q), (5) 


with 


V(p,q;M) = 


i as D„„(k) 7 f 72 + rl,(k)r?(k)r 2i y-k) + vl,(k) 


Here a s is the QCD coupling constant, D^ v is the gluon propagator in the Coulomb gauge 


4:71 
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F* (k) = — 777 , D'\k) = -- W- 
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k 2 ’ k 2 

and k = p — q; and u(p) are the Dirac matrices and spinors 

/ 


u x (p) = 


e(p) + m 

\ Hp) 


i 

crp 


X 


V e (p) + m 


( 6 ) 


(7) 


with e(p) = a/p 2 + m 2 . 

The effective long-range vector vertex of the quark is defined [16] by 

r M (k) =7^ +k = (0, k), (8) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In the configuration space the vector and scalar confining potentials in 
the nonrelativistic limit reduce to 


Vl,(r) = (1 - eJK^fr), 

KLM = eVUfM, (9) 

with 

Kont(r) = KLM + KUr) = At + B, (10) 

where e is the mixing coefficient. 


1 It is important to study diquark correlations in gauge-invariant color-singlet hadron states on the lattice. 
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(b) for diquark-antidiquark {dd!) interaction 

<d(P)|J M |d(Q)>4 


v (p> q; m) = 


otsD^i k) 


(d'{P')\J„\d'{Q')) 


2 V E d E d 3 2^/Ed* Ed 

+r d {P)rAP') \J*,„JSKA k) + KL(k)l MQ)MQ'), (11) 


where (d(P)\J li \d(Q)) is the vertex of the diquark-gluon interaction which takes into account 
the finite size of the diquark and is discussed below P ^ = (E d o , ±p) and = (E d o, ±q), 
E d = (M 2 - M 2 + MJ)/(2M) and E d , = (M 2 - M 2 + M 2 )/(2M) '. 

The diquark state in the confining part of the diquark-antidiquark quasipotential ED is 
described by the wave functions 


V’d(p) = 



for scalar diquark 

for axial vector diquark ’ 


where the four vector 

£d(p) = 


( (fjP) , (grfP)P 
V M d ’ d M d (E d (p ) + M d ) 


£d(p)P» = °> 


( 12 ) 


(13) 


is the polarization vector of the axial vector diquark with momentum p, E d (p ) = \Jp 2 + M d 
and e,i(0) = (0, s d ) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 



(P + Qh 


2y E d E d 
_ (.P + Q )„ 

2 J E d E d 


, JpAyul, 

+ 2JT7i 'I' ' 


for scalar diquark 
for axial vector diquark 


(14) 


where k = (0, k). Here the antisymmetric tensor 

(E^ = -i(^-^) (15) 

and the axial vector diquark spin is given by ( S d -k)u = — i£ku ; Pd is the total chromo- 
magnetic moment of the axial vector diquark. 

The constituent quark masses uif, = 4.88 GeV, m c = 1.55 GeV, m u = m d = 0.33 GeV, 
m s = 0.5 GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.3 
GeV have the usual values of quark models. The value of the mixing coefficient of vector 
and scalar confining potentials £ — — 1 has been determined from the consideration of 
charmonium radiative decays Hi and the heavy quark expansion Q. The universal Pauli 
interaction constant k — — 1 has been fixed from the analysis of the fine splitting of heavy 
quarkonia 3 Pj- states [li|. In this case the long-range chromomagnetic interaction of quarks 
vanishes in accord with the flux tube model. Finally, we choose the total chromomagnetic 
moment of the axial vector diquark p d = 0. 2 Such a choice appears to be natural, since the 


2 In our previous papers jl3, lla] on heavy baryons we adopted the value = 2 following the analogy with 
QED, which might be valid only for the perturbative one-gluon exchange part of the quark potential. 
However, in the above cases the long-range chromomagnetic diquark moment always contributed being 
multiplied by the vanishing total chromomagnetic long-range quark moment and, thus its value was 
irrelevant for the numerical results. 
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TABLE I: Masses M and form factor parameters of heavy-light diquarks. S and A denote scalar 
and axial vector diquarks antisymmetric [Q, q] and symmetric { Q , q} in flavour, respectively. 


Quark 

content 

Diquark 

type 


Q = c 



-0 

II 

Q> 


M (MeV) 

£ (GeV) 

C (GeV 2 ) 

M (MeV) 

£ (GeV) 

C (GeV 2 ) 

[Q,q] 

S 

f973 

2.55 

0.63 

5359 

6.10 

0.55 

{Q,q} 

A 

2036 

2.51 

0.45 

5381 

6.05 

0.35 

[Q,s] 

S 

2091 

2.15 

1.05 

5462 

5.70 

0.35 

{Q,s} 

A 

2158 

2.12 

0.99 

5482 

5.65 

0.27 


long-range chromomagnetic interaction of diquarks proportional to jid then also vanishes in 
accord with the flux tube model. 

At a first step, we calculate the masses and form factors of the heavy-light cliquarks. As 
it is well known, the light quarks are highly relativistic, which makes the v/c expansion 
inapplicable and thus, a completely relativistic treatment is required. To achieve this goal 
in describing heavy-light cliquarks, we closely follow our recent consideration of the spectra 
of light cliquarks in heavy baryons and adopt the same procedure to make the relativistic 

23 ). 


(see discussion in Ref. 


quark potential local by replacing E\p (p) = yjm 2 2 + P 2 — 1> Fi ,2 
The resulting light-quark-heavy-quark interaction potential is the same as the light quark- 
quark interaction [18] and is equal to 1/2 of the Qq interaction in the heavy-light meson. 3 
We solve numerically the quasipotential equation with this complete relativistic potential 
which depends on the diquark mass in a complicated highly nonlinear way. The obtained 
ground-state masses of scalar and axial vector heavy-light cliquarks are presented in Table HI 
In order to determine the cliquark interaction with the gluon held, which takes into ac¬ 
count the cliquark structure, it is necessary to calculate the corresponding matrix element 
of the quark current between cliquark states. This diagonal matrix element can be parame¬ 
terized by the following set of elastic form factors 
(a) scalar cliquark (S') 


(S(P)|J„|S(g)>=A + (*: 2 )(C > + Q)„. 


(16) 


(b) axial vector cliquark (A) 


(A(P)\JMm = -le* d (P)-e d (Q)]hi(k 2 )(P + Q)„ 

+M* 2 ) {K(P) • Q]e di „(Q) + M Q) ■ P]elJP)} 

+M * 2) Aff K( P ) ' Q][£d{Q) ' P|(P + Q) "’ 


(17) 


where k = P — Q and £d{P) is the polarization vector of the axial vector cliquark (THU) . 

Using the quasipotential approach with the impulse approximation for the vertex function 
of the quark-gluon interaction we find 0 


h + (k 2 ) = h (k 2 ) = h 2 {k 2 ) = F{ k 2 ), 
h(k 2 ) = 0, 


3 The masses of the ground state heavy-light mesons are well reproduced with this Qq potential. 
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FIG. 1: 


The form factors F(r) for {c, q} (solid line) and {b, q} (dashed line) axial vector diquarks. 


F(k 2 ) = 


VE d M d f d 3 p 

Ed + Md 


(2vr) 3 

pk 


P + 


2 e 2 (p) 


€i (p) + m i 


€i(p + k) + €i(p) 


’\/ e i(p + k)e 1 (p)(ei(p) + mi) 


E d + M d ) \ ei(p + fc) + mi |.2^ei(p + fc)ei(p) 

*d( P ) + (l<- 2 ), (18) 


where 40 are the diquark wave functions. We calculated the corresponding form factors 
F(r)/r which are the Fourier transforms of F(k 2 )/k 2 using the diquark wave functions 
found by numerical solving the quasipotential equation. Our estimates show that this form 
factor can be approximated with a high accuracy by the expression 

F(r) = l-e- 5r - ?r2 , (19) 

which agrees with previously used approximations E3- The values of parameters £ and ( 
for heavy-light scalar diquark [Q, q] and axial vector diquark { Q , q } ground states are given 
in Table |U In Fig. m we plot the functions F(r) for {Q, q} axial vector cliquarks. 

At a second step, we calculate the masses of heavy tetraquarks considered as the bound 
state of a heavy-light diquark and antidiquark. For the potential of the heavy diquark- 
antidiquark interaction m we get 


V(r) = Vcoui ( r ) + Konf(r) 


1 


E\ E 2 

F onl (r) + ^ (A + _ 

Vd(E x F 2 \C(r) T/C! q i 

tU T^WJ 

ly ou] (r) - VS.„,(r) + ^ E ' & - 


1 

"3 



^AKonf (r) + Vq oul (r) — 


l4EiE2_ /1 /V , V _ „v■ , A 

4 M x M 2 \r confl J confl ’) 


?(Sir)(S 2 r) 
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TABLE II: Masses of charm diquark-antidiquark states (in MeV). S and A denote scalar and axial 
vector diquarks. 


State 

Diquark 


Mass 


J PC 

content 

cqcq 

cscs 

cscq/cqcs 

IS 





0++ 

SS 

3812 

4051 

3922 

l +± 

(SA±SA)/V 2 

3871 

4113 

3982 

0++ 

AA 

3852 

4110 

3967 

1+- 

AA 

3890 

4143 

4004 

2++ 

AA 

3968 

4209 

4080 

IP 






SS 

4244 

4466 

4350 

TABLE III: Masses of bottom diquark-antidiquark states (in 

MeV). S and A denote scalar and 

axial vector diquarks. 




State 

Diquark 


Mass 


J PC 

content 

bqbq 

bsbs 

bsbq/bqbs 

IS 





0++ 

SS 

10471 

10662 

10572 

l +± 

{SA±SA)/V 2 

10492 

10682 

10593 

0++ 

AA 

10473 

10671 

10584 

1+- 

AA 

10494 

10686 

10599 

2++ 

AA 

10534 

10716 

10628 

IP 





1“ 

SS 

10807 

11002 

10907 


+ I[ 

M s,s,J, 


(20) 

where 

^coui(r) = 

4 F 1 (r)F 2 (r 

) 



Cl o 

3 r 



is the Coulomb-like one-gluon exchange potential which takes into account the finite sizes 

of the diquark and antidiquark through corresponding form factors Fi^{r) 

Here Sq2 and L 


are spin operators of the diaquark and antidiquark and their orbital momentum. 

The diquark-antidiquark model of heavy tetraquarks predicts in US the existence of a 
flavour SU(3) nonet of states with hidden charm or beauty (Q = c,b ): four tetraquarks 
{[Qq\[Qq], q = u, cl) with neither open or hidden strangeness, which have electric charges 
0 or ±1 and isospin 0 or 1; four tetraquarks ([<3s][Q<?] and [<5g][Qs], q = u,d) with open 
strangeness (S = ±1), which have electric charges 0 or ±1 and isospin one tetraquark 
([<5s][Qs]) with hidden strangeness and zero electric charge. Since in our model we ne¬ 
glect the mass difference of u and d quarks and electromagnetic interactions, corresponding 


















TABLE IV: Thresholds for open charm decays and nearby hidden-charm thresholds. 


Channel 

Threshold (MeV) 

Channel 

Threshold (MeV) 

Channel 

Threshold (MeV) 

D°D° 

3729.4 

DfD~ 

3936.2 

D°Df 

3832.9 

D + D~ 

3738.8 

r/'J/i/j 

4054.7 

D ± Df 

3837.7 

D°D*° 

3871.3 

DfDf F 

4080.0 

D*°Df 

3975.0 

pj/i) 

3872.7 


4116.4 

D°D* ± 

3976.7 

D ± D*+ 

3879.5 

D*+Dr 

4223.8 


3988.6 


3879.6 



K*°J/ip 

3993.0 

jj* 0 jj*0 

4013.6 



£)*0]j*± 

4118.8 


TABLE V: Thresholds for open bottom decays. 


Channel 

Threshold (MeV) 

Channel 

Threshold (MeV) 

Channel 

Threshold (MeV) 

BB 

10558 


10739 

BB S 

10649 

BB* 

10604 

BfB * s T 

10786 

B*B S 

10695 

B*B* 

10650 

B*+B* s ~ 

10833 

B*B* 

10742 


tetraquarks will be degenerate in mass. A more detailed analysis [11] predicts that such mass 
differences can be of a few MeV so that the isospin invariance is broken for the [Qq , ][Q?] 
mass eigenstates and thus in their strong decays. The (non)observation of such states will 
be a crucial test of the tetraquark model. 

The calculated heavy tetraquark masses are presented in Tables [H] and IIIII The corre¬ 
sponding open charm and bottom thresholds are given in Tables ITVl and El We find that all 
S -wave tetraquarks with hidden bottom lie considerably below open bottom thresholds and 
thus they should be narrow states which can be observed experimentally. This prediction 
significantly differs from the molecular picture where bound B — B* states are expected 
to lie very close (only few MeV below) to the corresponding thresholds. 

The situation in the hidden charm sector is considerably more complicated, since most of 
the tetraquark states are predicted to lie either above or only slightly below corresponding 
open charm thresholds. This difference is the consequence of the fact that the charm quark 
mass is substantially smaller than the bottom quark mass. As a result the binding energies 
in the charm sector are significantly smaller than those in the bottom sector. 

In Table IVTl we compare our results for the charm diquark-antidiquark bound states with 


the predictions of Ref. 


11 


The differences in some of the mass values can be attributed 


to the substantial distinctions in the used approaches. We describe the diquarks as quark- 


quark bound systems dynamically and calculate their masses, while in Ref. [llj they are 


treated only phenomenologically. Then we consider the diquark-antidiquark interaction and 
the tetraqark as purely the diquark-antidiquark composite system. In distinction Maini et 
ah consider a hyperhne interaction between all quarks which, e.g., causes the splitting of 
1 ++ and l + ~ states arising from the SA diquark-antidiquark compositions. From Table lYll 
where we also give possible experimental candidates for the neutral tetraquarks with hidden 


charm, we see that our calculation supports the assumption of Ref. [ll| that 77(3872) can 


be the axial vector 1 ++ tetraquark state composed from scalar and axial vector diquark 
and antidiquark in the relative S'-wave state. On the other hand, in our model the lightest 
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TABLE VI: Comparison of theoretical predictions for the masses of charm diquark-antidiquark 
states cqcq (in MeV) and possible experimental candidates. 


State 


Theory 


Experiment [^, 4, 

5, 6, 7, 21, 22, 23, 24] 

J PC 

this work 

□JJ 

[12] ( cscs ) 

state 

mass 

IS 






0 ++ 

3812 

3723 




1 ++ 

3871 

3872t 


A (3872) 

3871.2 ± 0.4 

1 +- 

3871 

3754 




0 ++ 

3852 

3832 




1 +- 

3890 

3882 




2 ++ 

3968 

3952 


Y (3943) 

3943 ±11 ±13 

IP 






1— 

4244 


4330 ± 70 

Y (4260) 

4259 ± 8± 2 6 


t input 


scalar 0 ++ tetraquark is predicted to be above the open charm threshold DD and thus to 
be broad, while in the model of Ref. |ll| it lies a few MeV below this threshold, and thus is 
predicted to be narrow. Our 2 ++ tetraquark also lies higher than the one in Ref. h|, thus 
making the interpretation of this state as Y (3943) less probable but still compatible with 
experiment. We also find that Y'(4260) cannot be interpreted as P -wave 1 state of charm- 
strange diquark-antidiquark, since its mass is found to be ~ 200 MeV heavier (see Table 111)) . 
A more natural tetraquark interpretation could be the P -wave ([cg]s = o[cg]s=o)p- W ave state 
which mass is predicted in our model to be close to the mass of V(4260) (see Table IYT1) . 
Then the dominant decay mode of V(4260) would be in DD pairs. 

In summary, we calculated the masses of heavy tetraquarks with hidden charm and bot¬ 
tom in the diquark-antidiquark picture. In contrast to previous phenomenological treatments 
we used the dynamical approach based on the relativistic quark model. Both diquark and 
tetraquark masses were obtained by numerical solution of the quasipotential equation with 
the corresponding relativistic potentials. The diquark size was also taken into account with 
the help of the diquark-gluon form factor in terms of diquark wave functions. It is important 
to emphasize that, in our analysis, we did not introduce any free adjustable parameters but 
used their fixed values from our previous considerations of heavy and light meson proper¬ 
ties. It was found that the A"(3872) can be the neutral charm tetraquark state. If it is 
really a tetraquark, one more neutral and two charged tetraquark states must exist with 
close masses. The ground states of bottom tetraquarks are predicted to have masses below 
the open bottom threshold and thus should be narrow. The (non)observation of these states 
will be an important test of the tetraquark model. 
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